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Abstract. We present evolutionary calculations for 74 close binaries systems with initial primary masses in the range 
12.. .25 M0, and initial secondary masses between 6 and 24 M0. The initial periods were chosen such that mass overflow starts 
during the core hydrogen burning phase of the primary (Case A), or shortly thereafter (Case B). We use a newly developed 
binary code with up-to-date physics input. Of particular relevance is the use of OPAL opacities, and the time-dependent 
treatment of semiconvective and thermohaline mixing. We assume conservative evolution for contact-free systems, i.e., no 
mass or angular momentum loss from those system except due to stellar winds. 

We investigate the borderline between contact-free evolution and contact, as a function of the initial system parameters. The 
fraction of the parameter space where binaries may evolve while avoiding contact — which we found already small for the 
least massive systems considered — becomes even smaller for larger initial primary masses. At the upper end of the considered 
mass range, no contact-free Case B systems exist. While for primary masses of 16 Mq and higher the Case A systems dominate 
the contact-free range, at primary masses of 12 M0 contact-free systems are more frequent for Case B. We identify the drop of 
the exponent x in the main sequence mass-luminosity relation of the form L oc as the main cause for this behaviour 
For systems which evolve into contact, we find that this can occur for distinctively different reasons. While Case A systems 
are prone to contact due to reverse mass transfer during or after the primary's main sequence phase, all systems obtain contact 
for initial mass ratios below ~ 0.65, with a merger as the likely outcome. We also investigate the effect of the treatment of 
convection, and found it relevant for contact and supernova order in Case A systems, particularly for the highest considered 
masses. 

For Case B systems we find contact for initial periods above ^ 10 d. However, in that case (and for not too large periods) 
contact occurs only after the mass ratio has been reversed, due to the increased fraction of the donor's convective envelope. 
As most of the mass transfer occurs conservatively before contact is established, this delayed contact is estimated to yield to 
the ejection of only a fraction of the donor star's envelope. Our models yield the value of /3, i.e., the fraction of the primaries 
envelope which is accreted by the secondary. 

We derive the observable properties of our systems after the major mass transfer event, where the mass gainer is a main 
sequence or supergiant O or early B type star, and the mass loser is a helium star We point out that the assumption of 
conservative evolution for contact-free systems could be tested by finding helium star companions to O stars. Those are also 
predicted by non-conservative models, but with different periods and mass ratios. We describe strategies for increasing the 
probability to find helium star companions in observational search programs. 
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1. Introduction 

To understand the evolution of massive close binary system is 
relevant for many current astrophysical problems. Massive bi- 
naries provide the unique means to obtain accurate fundamen- 
tal stellar properties of massive stars for various evolutionary 
stages. I.e., they have been used to derive masses of Wolf-Rayet 
stars (van der Hucht 2000) and to establish mass-luminosity 
relations for massive stars (Martin et al 1998, Ostrov et al. 
2000). They also lead to the existence of stellar configurations 
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which do not evolve from single stars and which have highly 
exciting observational manifestations. Examples are the pecu- 
liar supernova 1987 A, which several authors attribute to binary 
evolution (Podsiadlowski 1992, de Loore & Vanbeveren 1992, 
Braun & Langer 1995), short-period Wolf-Rayet binaries (e.g., 
Harries & Hilditch 1997), massive X-ray binaries (Chevalier & 
Ilovaisky 1998), supernovae of type lb and Ic (Podsiadlowski 
et al. 1992), and 7-ray burster (Fryer et al. 1999) 

However, the evolutionary theory of massive close bina- 
ries still suffers from fundamental uncertainties. The evolution 
of such systems has been modelled by various authors, e.g., 
Paczynski (1967), Kippenhahn (1969), de Loore & De Greve 
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(1992), Pols (1994), Vanbeveren et al. (1998abc), Wellstein & 
Langer (1999). One longstanding important question in binary 
models which include mass transfer is: how much of the trans- 
fered matter can actually be accreted by the secondary star? In 
most calculations so far, the fraction (3 of the accreted matter 
which is retained on the accreting star is chosen as a constant, 
with P = 0.5 being a frequent choice (e.g., Moers & van den 
Heuvel 1989, De Greve & De Loore 1992, De Loore & De 
Greve 1992, Vanbeveren et al. 1998abc). This means that half 
of the matter lost by the donor star is accreted by the mass 
gainer; the other half supposed to leave the binary system due 
to an undefined force, and carrying an amount of angular mo- 
mentum which needs to be described by a second parameter 
(cf. Huang & Taam 1990, Podsiadlowski et al. 1992). Also the 
assumption of a fully conservative mass transfer where all mat- 
ter is accreted by the secondary, i.e. (3—1, has been used. 
Physical models which yield the parameter /? as a function of 
the system parameters, or even as a function of time for a give 
system, are lacking. 

Arguments of /3 ~ 0.5 on average are being used in the 
literature (De Greve & de Loore 1992, de Loore & De Greve 
1992, Vanbeveren et al. 1998abc), and correspondingly mas- 
sive close binary models using /3 = 0.5 have been computed 
(see above). Here, we pursue a different approach. We adopt 
/? = 1 as long as the two stars in the binary system do not 
evolve into contact. We want to find out how many massive bi- 
naries actually avoid contact when /? = 1 is used, and whether 
the result can be compatible with an average value of /3 ~ 0.5. 

Mass transfer in systems with initial mass ratios very differ- 
ent from unity is known to lead to the rapid expansion of the ac- 
creting star and thus to a contact system (Benson 1970, Ulrich 
& Burger 1976). Also systems with large initial periods are 
prone to develop contact, as convective envelopes of the donor 
stars lead to very high mass transfer rates (cf. Podsiadlowski et 
al. 1992). Pols (1994) found that Case A[| systems with primary 
masses in the range 8... 16 M0 avoid contact for initial mass ra- 
tios g ^ 0.7 anc/ initial periods Pq > 1.6d. However, Pols used 
opacities from Cox & Stewart (1970), the Schwarzschild crite- 
rion for convection, and instantaneous thermohaline mixing. In 
the present work, we reconsider the work done by Pols using 
updated physics input, and we extend the considered parame- 
ters space towards higher masses and larger initial periods. 

We describe our computational method in Section 2, and 
illustrate our binary model computation at the example of two 
representative cases in Section 3. In Section 4, we present the 
results of 74 model systems and discuss the formation of con- 
tact as a function of the main system parameters and the applied 
stellar physics. We discuss our results in Section 5, attempt a 
comparison with observations in Section 6, and summarise our 
main conclusions in Section 7. 



' Case A: mass transfer during central hydrogen burning; Case B: 
mass transfer after the end of central hydrogen burning; Case C: mass 
transfer after the end of central helium burning 



2. Computational methods 

We computed the evolution of massive close binary systems us- 
ing a computer code generated by Braun (1997) on the basis of 
an implicit hydrodynamic stellar evolution code for single stars 
(cf. Langer 1991, 1998). It invokes the simultaneous evolution 
of the two stellar components of a binary and computes mass 
transfer within the Roche approximation (Kopal 1978). 

Mass loss from the Roche lobe filling component through 
the first Lagrangian point is computed according to Ritter 
(1988) as 

^ 

M = Mo exp with Mq = ^pVsQ. (1) 

Hp 

Hp is the photospheric pressure scale height, p the density, 
Vs the velocity of sound and Q the effective cross-section of 
the stream through the first Lagrange pint according to Meyer 
& Meyer-Hofmeister (1983). The distance between the first 
Lagrange point and the center of mass which which enter 
Hp and Q are interpolated from the tables of Kopal (1978) 
and Mochnacki (1984). We note that this treatment is only 
correct for stars with small photospheric pressure scale com- 
pared to their radius (Pastetter & Ritter 1989). Otherwise, 
models of Paczyriski and Sienkiewicz (1972), Plavec et al. 
(1973), Savonije (1978) orKolb & Ritter (1990) might be used. 
However, this is not required for the binaries considered in this 
work. 

Due to the exponential dependence of the mass transfer 
rate on the stellar radius, an explicit scheme for computing 
the mass transfer rate can easily lead to numerical instabilities. 
Thus, we derive the mass transfer rate implicitly using an it- 
erative scheme (combined Secant/Bisection-method; cf. Press 
et al. 1988), which is numerically more stable (Braun 1997). 
However, for each iteration a complete stellar model needs to 
be computed. 

The change of the orbital period due to the mass trans- 
fer and stellar wind mass loss is computed according to 
Podsiadlowski et al. (1992), with the specific angular momen- 
tum of the stellar wind material being determined according to 
Brookshaw & Tavani (1993). The spin angular momentum of 
both stars is neglected. Mass transfer is treated quasi conserva- 
tively. The only mass loss from the system is due to the stellar 
winds of both components. The adopted stellar wind mass loss 
rates are described in Wellstein & Langer (1999). 

We include time-dependent thermohaline mixing — i.e. 
mixing in radiatively stable regions which occurs due to an 
outwards increase of the mean molecular weight. This mix- 
ing process is important in stars which accrete helium-enriched 
material. A time-dependent treatment of thermohaline mixing 
is important for identifying the borderline between contact and 
contact-free close binary systems as it occurs on the thermal 
time scale of the accreting star, which — at the borderline of 
contact-free systems — is just equal to the mass transfer time 
scale. I.e., a time dependent treatment is required to reproduce 
the radius evolution of the accreting star during the mass trans- 
fer phase correctly as soon as helium-enriched material is being 
transferred (cf., Braun 1997). 
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Numerically, we treat thermohaline mixing through a dif- 
fusion scheme (Braun 1997). The corresponding diffusion co- 
efficient is based on the work of Stern (1960), Ulrich (1972), 
and Kippenhahn et al. (1980); it reads 



D 



thn 
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2cpp \/ad - V 

where K — 4ac/{3Kp), (j> — {d\np/d\nfi)pT, 

5 = - (a In p/a In T)p,^, = dln^l/dlnP, V ad = 
[d\ivT/d\nP)ad, and V = dlnT/dlnP. The quantity at,„„ 
is a efficiency parameter of order unity, and we use athm = 1 
throughout this paper. 

The entropy of the accreted material is assumed to be equal 
to that of the surface of the mass gainer and the gravitational 
energy release due to mass transfer is treated as in Neo et 
al. (1977). Convection and semiconvection are treated as de- 
scribed in Langer (1991) and Braun & Langer (1995) (c.f. also 
Langer et al. 1983). A semiconvective efficiency parameter of 
Use = 0.01 is applied in most models. We do not include so 
called convective core overshooting in our models, as recent 
evidence impUes that the physical effect of increased convec- 
tive core masses in massive main sequence stars may be due to 
effects of rapid rotation rather than convection (Maeder 1987, 
Langer 1992, Heger et al. 2000). As the rotation rates in close 
binaries may be reduced due to orbit circularisation and spin- 
orbit synchronisation (e.g., Savonije & Papaloizou 1997), it 
is unclear to what extent any convective core mass increase 
would occur in those objects. Opacities are taken from Iglesias 

6 Rogers (1996). Changes in the chemical composition are 
computed using a nuclear network including the pp-chains, the 
CNO-tri-cycle, and the major helium-, carbon, neon and oxy- 
gen burning reactions. Further details about the computer pro- 
gram and input physics can be found in Langer (1998) and 
Wellstein & Langer (1999). For all models, a metallicity of 2% 
is adopted. 

3. Case A and B evolution: examples 

We have computed 74 Case A and Case B systems to derive 
the limits of contact free evolution for massive close binaries 
within the assumptions made in Section 2. Before we present 
the results of these models, we want to discuss the evolution of 
one typical Case A system (No. 3 1 ; cf. Sect. 4 below), and one 
Case B system (No. 37) in the following, to illustrate the capa- 
bilities of our method and the type and quality of the derived re- 
sults. Those two systems are picked to show one detailed exam- 
ple per considered initial primary star mass — i.e., 12 Mq and 
16 Mq; we have discussed details of models with 25 ini- 
tial primary mass in Wellstein & Langer (1999). Furthermore, 
the two examples shown in detail below correspond to systems 
which evolve contact-free. Therefore, we could follow the evo- 
lution of both binary components until the end of carbon burn- 
ing (where we stopped the calculations). 



3.1. Case A 

Table [l| and Fig. |l] describe the evolution of our Case A sys- 
tem No. 31 (c.f. Table |] below), which initially consists of a 



12 M0 primary ^ and a 7.5 Mq secondary star in a 2.5 d orbit. 
Mass transfer starts at a system age of ^ 1.3 10^ yr, at a core 
helium mass fraction of the primary of 0.89. The mass trans- 
fer rate rises to some 10^"* Mq yr^^ during the so called rapid 
Case A mass transfer phase (Fig. i.e. to values of the or- 
der of A/i/tkh,!, Ml being the primary mass and tkh,! the 
thermal time scale of the primary. The reason is that the orbital 
separation of both stars decreases until both components have 
the same mass, here about 9.75 Mq. During the slow Case A 
mass transfer phase, the mass transfer is driven by the nuclear 
evolution of the primary, and the mass transfer rate drops to 
some 10^^ Mq yr^^, i.e. values of the order of Afi/r„uc,i- 

It is remarkable that the mass transfer does not immedi- 
ately relax at the time of minimum orbital separation. In fact. 
Fig. ^ shows that only a small amount of mass 0.15 Mq) 
is transfered during that slow Case A mass transfer, i.e. that 
its rapid phase does not stop until the primary mass is as 
small as 4Mq. The reason is that, in contrast to the rule of 
thumb that stars with radiative envelopes shrink upon mass 
loss, the evolved core hydrogen burning 12 Mq star has a neg- 
ative mass-radius exponent, i.e. its thermal equilibrium radius 
increases for decreasing mass as a response to the increasing 
core-mass/envelope-mass ratio. The thermally unstable mass 
transfer during the rapid Case A is ended only when the mass 
ratio becomes sufficiently small (cf. Ritter 1988). 

During the rapid phase of the Case A mass transfer the 
primary becomes underluminous and the secondary becomes 
overluminous, due to the rapid mass loss and gain, respectively. 
I.e., they deviate by up to a factor of two from their thermal 
equilibrium luminosity (see the kinks in the evolutionary tracks 
towards the end of the Case A mass transfer in Fig. |l]). During 
the slow phase of the Case A mass transfer, primary and sec- 
ondary relax to thermal equilibrium. Due to the change in total 
mass, the secondary becomes more luminous than the primary, 
even though the primary, after the mass transfer, is strongly 
overluminous, i.e. its luminosity exceeds that of a 4 Mq single 
star by a factor of 10. 

At core hydrogen exhaustion, the primary mass is about 
3.9 Mq. As it still contains a hydrogen-rich envelope, a hy- 
drogen burning shell ignites and the star attempts to expand 
to red giant dimensions. This initiates the so called Case AB 
mass transfer, which proceeds on the thermal time scale of 
the primary and produces correspondingly high mass transfer 
rates (cf. Fig. 2). It leads to the loss of almost all remaining 
hydrogen and to the extinction of the hydrogen burning shell. 
The Case AB mass transfer reduces the primary mass to about 
1.45 Mq, which coiTesponds to the mass of the convective hy- 
drogen burning core at central hydrogen exhaustion. Due to 
the extreme mass ratio (cf. Table 1) the orbital period grows 
from 8 d to about 100 d during the Case AB mass transfer. The 
primary spends about 5 Myr for core helium burning, during 
which it remains on the helium-main sequence in the HR dia- 
gram (cf. Fig. 1). The ensuing helium shell burning expands the 



We designate the initially more massive star in the system as the 
primary, also during its later evolution when it has become the less 
massive and/or less luminous star. The other star is designated as the 
secondary. 



4 



S. Wellstein et al.: Formation of contact in massive close binaries 




4.8 4.6 4.4 4.2 4.0 

log(T.„/K) 

Fig. 1. Evolutionary tracks of the primary (solid and dotted Une) and secondary star (dashed line) of our case A binary system 
No. 31 (initial masses are 12 Mq and 7.5 M0, the initial period is 2.5 d; of. Table 3 in Section 4) in the HR diagram. Beginning 
and the end of the mass transfer phases are marked with numbers; 1: begin of Case A, 2: end of Case A, 3: begin of Case AB, 
4: end of Case AB. The labels A/a designate the end of central hydrogen burning of the primary/secondary, B/b the end of 
central helium burning of the primary/secondary, and c the point of the supernova explosion of the secondary. In this system, 
the secondary star ends its evolution first. The time of its supernova explosion marks the end of the solid line in the track of the 
primary. The further evolution of the primary is shown as dotted Une. During this phase, it is treated as a single star since the 
system is likely broken up due to the secondary's explosion. 




time [yr] 

Fig. 2. Mass transfer rate as function of time for our Case A 
system No. 3 1 (initial masses are 12 Mq and 7.5 Mq, the initial 
period is 2.5 d; cf. Table 3 in Section 4) The left panel shows 
the rapid and slow Case A mass transfer, the right panel shows 
the ensuing Case AB mass transfer. 



star to red giant dimensions, and a so called Case ABB mass 
transfer would occur (Kippenhahn & Thomas 1979) — it does 
occur in other systems of our sample — if the secondary star 
would not have exploded in the meantime (see below). We fol- 
low the primary's evolution until core carbon exhaustion. Due 
to its small final mass, it is unlikely to produce a supernova 
explosion. Instead, mass loss during its red giant stage at mod- 
erate rates would suffice to substantially reduce the total mass 
before the core mass can grow to ignite neon burning (Habets 
1986ab, Woosley et al. 1995). I.e., our primary star of initially 
12 Mq ends as an isolated ONeMg- white dwarf. 

The secondary star grows to 17.81 Mg, and finishes core 
hydrogen burning about 3.5 Myr after the end of the case AB 
mass transfer. At this time it begins to expand to become a su- 
pergiant. However, due to the presence of a strong mean molec- 
ular weight barrier above the convective core at the onset of 
the mass transfer, it does not rejuvenate during core hydro- 
gen burning, i.e. it does not adapt its convective core mass to 
its new total mass (cf. Hellings 1983, 1984; Braun & Langer 
1995). Therefore, the secondary star has, after the mass trans- 
fer, a much smaller helium core to envelope mass ratio as a 
single star of the same mass. For this reason, it remains in the 
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Fig. 3. Evolutionary tracks of the primary (solid and dotted Une) and secondary star (dashed line) of our case B binary system 
No. 37 (initial masses are 16 M0 and 15 Mq, the initial period is 8 d; of. Table 3 in Section 4) in the HR diagram. Beginning and 
end of the mass transfer phases are marked with numbers; 1: beginning of case B, 2: end of case B, 3: beginning of case BB, 4: 
end of case BB. The labels A/a designate the end of central hydrogen burning of the primary/secondary, B/b the end of central 
heUum burning of the primary/secondary, and C/c the point of their supernova explosion. 



Table 1. Evolutionary parameters of the Case A system No. 3 1 . 
The calculation ends after core carbon exhaustion of the pri- 
mary (note that the secondary star produces the first super- 
nova in the system). The different columns have the following 
meaning: evolutionary stage (ZAMS = zero age main sequence; 
ECHB = end core hydrogen burning, ECHeB = end core he- 
lium burning; ECCB = end core carbon burning), evolutionary 
age (ZAMS: t = 0), primary and secondary star mass, orbital 
period, primary and secondary star orbital velocity. 





time 


Ml 


M2 


P 


Vl 


V2 




M3T 


Mo 


Mo 


d 


km/s 


km/s 


ZAMS 





12.0 


7.50 


2.5 


162 


260 


begin Case A 


12.96 


11.9 


7.49 


2.52 


162 


258 


end Case A 


14.93 


3.92 


15.5 


7.99 


228 


57.8 


ECHB primary 


15.17 


3.92 


15.5 


7.99 


228 


57.8 


begin Case AB 


15.17 


3.92 


15.5 


7.99 


228 


57.8 


end Case AB 


15.46 


1.45 


17.8 


102.6 


113 


9.2 


ECHB secondary 


18.94 


1.45 


17.7 


104.2 


112 


9.2 


ECHeB primary 


20.55 


1.45 


174 


108.0 


110 


9.2 


ECHeB secondary 


20.82 


1.45 


17.2 


108.1 


109 


9.2 


ECCB secondary 


21.16 


1.45 


17.0 


1084 


109 


9.2 


ECCB primary 


22.14 


1.45 











blue part of the HR diagram and does never become a red su- 
pergiant (Braun & Langer 1995). 

Whether or not a secondary star rejuvenates depends 
strongly on the — still poorly known — efficiency of semi- 
convection mixing (Langer et al. 1983, Braun & Langer 1995). 
Were rejuvenation occurring in our secondary, it would, after 
the mass transfer, evolve similar to a single star of ^ 18 M0 
and attempt to become a red supergiant after core hydrogen ex- 
haustion. This would lead to a reverse Case B mass transfer in 
this system, with the consequence of contact evolution and a 
hkely merging of both stars, due to the extreme mass ratio (cf . 
System No. 25 below, and Section 5.2, for more details). 

Because the secondary star has a small helium core mass, 
its fuel supply for core helium burning is small. However, due 
to its comparatively large total mass, its luminosity is large. 
Both factors together imply a very short core helium burning 
life time of only 2 Myr. As a consequence, the secondary over- 
takes the primary in the evolutionary status and ends its evo- 
lution first. It explodes as a nitrogen-rich blue supergiant with 
properties similar to those of the progenitor of SN 1987A. It 
is helium- and nitrogen enriched due to the accretion, by fac- 
tors 1.26 and 5.4, respectively, while carbon and oxygen are 
depleted by factors of 0.23 and 0.81 (cf. Wellstein & Langer 
2001). Its core still has a memory of the lower initial mass: the 
helium core mass is only ^ 4 Mq at core carbon exhaustion. 
This model applied to SN 1987A would predict the existence 
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time [yr] 

Fig. 4. Mass transfer rate as function of time for our Case B 
system No. 37 (initial masses are 16 M0 and 15 M0, the initial 
period is 8 d; cf. Table 3 in Section 4) The left panel shows 
the Case B mass transfer, the right panel shows the ensuing 
Case BB mass transfer. 

of a 1.45 M0 helium star with a luminosity of ^ 10"^ Lq and a 
temperature of ^ 55 000 K in the supernova remnant. 

In many of our Case A systems — always for those with a 
12 Mq primary — the secondary is the first star which reaches 
the supernova stage. This supernova reversal is discussed in 
more detail in Section 5.2. 

3.2. Case B 

Table I and Fig. | describe the evolution of our case B sys- 
tem No. 37 (cf. Table |), starting with a 16 primary and 
a 15 Mq secondary in an 8 d orbit. As the primary expands 
on a thermal time scale after core hydrogen exhaustion, the 
Case B mass transfer occurs on the thermal time scale of the 
primary, with mass transfer rates of the order of Mi/tkh.i (cf. 
Sect. 3.1). Both stars are far from thermal equilibrium during 
the whole mass transfer phase. However, also in this case, a fast 
and a somewhat slower Case B mass transfer phase can be dis- 
tinguished (Fig. Q), since the orbit shrinks at first and expands 
after the mass ratio is reversed. 

In the considered system, the mass transfer achieves about 
2 10^"^ M0 yr~^ during the Case B mass overflow. This im- 
plies an accretion time scale of the secondary of the same or- 
der as its Kelvin-Helmholtz time scale. Consequently, the sec- 
ondary swells considerably during the Case B mass transfer 
While here, contact is not quite achieved, for an initial period 
of 9 d instead of 8 d the primary's radius at the onset of mass 
transfer is slightly larger, the mass transfer rate as well, and the 
system does evolve into contact at this stage (cf. system No. 39 
in Sect. 4 below). 

After the Case B mass transfer, very little hydrogen is left 
in the primary star. It evolves quickly into a hot and compact 
helium star of ~ 3.8 Mq, as which it burns helium in its core 
for about 1 Myr. Thereafter, it expands again (cf. Sect. 3.1) 
and fills its Roche lobe a second time, which initiates the so 
called case BB mass transfer During this phase, it loses the 



Table 2. Evolutionary parameters of the Case B system No. 37. 
Both stars are followed beyond core carbon exhaustion. The 
different columns have the following meaning: evolutionary 
stage (ZAMS = zero age main sequence; ECHB = end core 
hydrogen burning, ECHeB = end core helium burning; ECCB 
= end core carbon burning), evolutionary age (ZAMS: t = 0), 
primary and secondary star mass, orbital period, primary and 
secondary star orbital velocity. 





time 


Ml 


Ma 


P 


^"1 


V2 




Myr 


M© 


Mo 


d 


km/s 


km/s 


ZAMS 





16.0 


15.0 


8 


162 


173 


ECHB primary 


9.472 


15.8 


14.8 


8.19 


160 


170 


begin Case B 


9.490 


15.8 


14.8 


8.19 


160 


170 


end Case B 


9.513 


3.83 


26.7 


94.9 


128 


18.3 


ECHeB primary 


10.46 


3.61 


26.5 


99.1 


126 


17.1 


begin Case BB 


10.68 


2.81 


26.5 


102.9 


126 


13.4 


ECCB primary 


10.68 


2.33 


26.8 


170.7 


109 


9.4 


end Case BB 


10.68 


2.32 


26.8 


172.1 


108 


9.4 


ECHB secondary 


10.95 




26.7 








ECHeB secondary 


11.71 




26.0 








ECCB secondary 


11.83 




25.9 









major fraction of its helium envelope — 0.49 Mq — to the 
secondary. This mass transfer occurs again on the thermal time 
scale of the primary (cf. Fig. Between core carbon exhaus- 
tion and neon ignition, the mass transfer ends. During the fur- 
ther evolution the binary system remains detached, but the pri- 
mary remains an extended supergiant and nearly fills its Roche 
lobe. In the final state, the once 16 Mq primary is a 2.3 Mq 
helium giant with a remaining helium mass of 0.74 Mq. The 
core of the primary will collapse and produce a supernova ex- 
plosion. The explosion characteristics may resemble those of 
Type lb or Ic supernovae (Shigeyama et al. 1990, Woosley et 
al. 1995). However, the helium envelope of our model star is 
very extended; its radius is ~ 80 Rq. Theoretical models for 
the explosions for such cool helium giants are still missing in 
the literature. Due to the large final mass ratio of 27/2.3, this 
system has a good chance to remain bound after the explosion 
of the primary component. 

The secondary star grows to almost 27 Mq during the 
Case B mass transfer. It continues its main sequence evolution 
thereafter and is still burning hydrogen in its core when the pri- 
mary explodes. During the Case BB mass transfer, it accretes 
0.49 Mq of almost pure helium from the primary. During this 
phase, it becomes very helium and nitrogen rich at the surface. 
Even though during and after the Case BB mass transfer this 
enrichment is strongly diluted by thermohaline mixing, the fi- 
nal helium and nitrogen enrichment factors are 1.25 and 4, i.e. 
the final surface mass fractions are 0.35 and 0.0041 for helium 
and nitrogen, respectively. 

After core hydrogen exhaustion, the secondary expands, 
but it remains a blue supergiant up to its supernova explosion 
(Fig. 3) due to the lack of rejuvenation (cf. Sect. 3.1). The ex- 
plosion of the secondary disrupts the binary system with a high 
probability. 
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4. Formation of contact 

In this section, we present results for 74 computed systems. 
Rather than providing a detailed description of the evolution of 
each system, we focus on the general trends of key properties 
as function of the major initial parameters. Table 3 provides 
the most important quantities for all computed systems. We fo- 
cus on the various mechanisms which can drive massive close 
binary systems into contact. 

4. 1. Contact formation due to small or large periods 
4.1.1. Case B 

Case B systems with an initial mass ratio close to 1 and a 
rather short initial period (early Case B) are known to have 
a good chance to avoid contact during the Case B mass transfer 
(cf. Pols 1994). Our Table 3 contains 15 examples, i.e. sys- 
tems No. 2. ..6, 12, 16, 18. ..20, 27, 28 (initial primary mass 
of 12 Mq), and No. 35, 37 and 44 (initial primary mass of 
16 Mq). The fact that no such system with an initial primary 
mass of 25 Mq could be found is discussed in Sect. 5.1. It is 
often assumed that late Case B systems, i.e. such with a rela- 
tively large initial period, evolve through a common envelope 
stage (Podsiadlowski et al. 1992, Vanbeveren 1998abc). In the 
following, we will have a closer look to the transition region in 
between both extremes. 

Fig. H shows the mass transfer rates as a function of the 
secondary mass for Case B binaries with initial primary and 
secondary masses of 12 Mq and 1 1 Mq, respectively, but dif- 
ferent initial periods in the range 6.. .40 d (Systems No. 3 to 8). 
Generally, initially wider systems develop larger mass trans- 
fer rates. However, it is interesting to consider the time depen- 
dence of the mass transfer rate of these systems in some detail. 
The first maximum in the mass transfer rate is related to the 
end of decrease of the orbital separation at mass ratio one. I.e., 
when the primary has become the less massive component in 
the system, the period and the Roche lobe of the primary be- 
gin to grow, and the increase of the mass transfer rate flattens. 
This happens with some time delay due to the finite thermal 
response time of the primary. 

Later during the mass transfer, the orbit widens signifi- 
cantly. Therefore, the primary becomes more extended and its 
surface temperature decreases. In the wider systems considered 
in Fig. ^ the outer envelope of the primary becomes even con- 
vectively unstable. It is the drop of the adiabatic mass-radius 
exponent (^ad (cf. Ritter 1988) which leads to another rise of 
the mass transfer rate even though the orbit widens. We desig- 
nate this feature as delayed contact. Table 3 shows nine systems 
which evolve delayed contact. It is restricted to the Case B and 
has no counterpart in our Case A systems (cf. Sect.4.1.2). 

Fig. ^ displays the radius evolution of the secondary stars 
of the systems considered in Fig. |[ Initially, the radii do not 
grow significantly, although it can be seen seen that larger 
radii are obtained for larger initial periods, i.e. for larger mass 
transfer rates. Later on, mass transfer rates in excess of ^ 
3 10^^ Mq yr^^ make the secondary swell enormously. This 
mass transfer rate corresponds roughly to M — RL/{GM), 
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Fig. 5. Mass transfer rates as a function of secondary mass for 
Case B systems initially consisting of a 12 Mq and all Mq 
star with initial periods in the range 6...40d (No. 3,4,5,6,7 and 
8). The systems with 30 and 40 days initial period form contact, 
at which time we stop the calculation. 
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Fig. 6. Radius of the secondary during mass transfer as a func- 
tion of its mass, for the same systems as in Fig. ||. In the systems 
with 30 and 40 days initial period the large increase in radius 
leads to contact, at which time we stop the calculation. 



i.e. to the condition M/M = GM'^/{RL) for the secondary. 
In the two systems with initial periods of 30 and 40 days, the 
radius increases until the secondary also fills its Roche lobe 
and a contact system is formed. This effect marks the limit of 
contact free systems towards larger periods. 

It is remarkable that approximately 5 Mq for the 40 d bi- 
nary (system No. 8) and 6 Mq for the 30 d binary (system 
No. 7) are transfered before the contact occurs (cf. Fig. 6). 
From Fig. || we see that contact occurs the earlier (in terms 
of already transfered mass) the larger the initial system period 
is. I.e., there is a continuous transition towards very wide bina- 
ries with convective primaries, which evolve into contact at the 
very beginning of the mass transfer process (Podsiadlowski et 
al. 1992). 
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Table 3. Characteristic quantities for all computed systems: system number, initial values of primary mass Mi, secondary mass 
M2, and period Pj, and the mass transfer case. Bold values for the initial masses identify the component which ends its evolution 
first. Further columns are: the kind of contact experienced by the system (no contact='-', premature contact='M', reverse con- 
tact='R', q-contact='Q', and delayed contact='D'; a lower case letter is used if the secondary fills less than 1.5 times its Roche 
lobe; systems No. 52 and 59 experience a second contact after a weak q-contact), the core hehum mass fraction of the primary at 
the onset of the first mass transfer Yi, and M[, M^, L[, L'2, R'l, R2, P' denote primary and secondary mass, luminosity, effective 
temperature, radius and orbital period after the Case A/AB or Case B mass transfer, at a time when the core helium mass fraction 
of the primary has decreased to 0.8 due to heUum burning. Mmax is the maximum mass transfer rate (lower limits are give for 
contact systems, values in italic are uncertain) and AM is the amount of mass transfered by the time the maximum mass transfer 
rate or contact is reached. designates the fraction of the transfered matter that is accreted, assuming no further accretion after 
the secondary overfills its Roche lobe more than 1 .5 times. 
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= Secondary at ECHeB 
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According to the estimate of Webbink (1984), which com- 
pares the orbital energy of the binary with the binding energy 
of the envelope, our System No. 7 does not merge during the 
Case B mass transfer. Instead, with an efficiency parameter of 
Q!cc = 1, it obtains a period of about 7.5 d in a common en- 
velope evolution by expelhng the remaining envelope of the 
primary of ^ 3 from the system. 

In summary, in delayed contact systems a conservative part 
of Case B mass transfer is followed by a non-conservative one, 
leading to an oscillation of the orbit: it first shrinks until the 
mass ratio is 1, then widens up to the end of the conservative 
part, and then shrinks again due to mass and angular momen- 
tum loss from the system. For our System No. 7, the parame- 
ter /? for the complete Case B mass transfer, i.e., the ratio of the 
amount of mass accreted by the secondary to the amount lost 
by the primary, is about /3 ~ 6 /9 Mq ~ 0.67. 

After the delayed contact and coimnon envelope evolution. 

System No. 7 consists of a ^ 20 Mq main sequence star — the 
secondary — and a ~ 2.4 Mq hehum star in a 7.5 d orbit. The 



further evolution (which we have not computed) depends on 
the post-common envelope period. For 7.5 d, the primary fills 
its Roche lobe as it expands to a hehum giant, and transfers 
part of its hehum envelope to the secondary in a Case BB mass 
transfer, similar to that discussed in Sect. 3.2 for system No. 37. 

In contrast, the systems which avoid contact during the 
Case B mass transfer all together can often, even if the sec- 
ondary evolves faster, avoid the reverse Case B mass transfer 
if the secondary star does not rejuvenate (cf. values of T^'jj 
in Table 3). After the Case B mass transfer, their periods are 
typically larger than 50 d. 

We find that in our most massive systems, with a primary 
mass of 25 Mq, even early Case B mass transfer leads into con- 
tact, even for an initial mass ratios very close to 1 (cf. sys- 
tems No. 68 and 69 in Table 3). In the two considered sys- 
tems, large amounts of mass are transferred conservatively be- 
fore contact occurs (8.2 Mq and 7.4 Mq; cf. Table 3), i.e., we 
have a delayed contact situation. In system No. 68, the contact 
is marginal, in system No. 69 it is not. 
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Fig. 7. Secondary radius as a function of its mass during Case A 
mass transfer, for three systems starting out with a 16 M0 and 
a 12 M0 component, with initial periods in the range 1.5. ..2 d 
(Systems No. 53. ..55). In the two closer systems, the secon- 
daries terminate core hydrogen burning during the slow phase 
of the Case A mass transfer Their ensuing rapid expansion 
leads to contact. In the 2 d period system No. 55 (solid line). 
Case AB mass transfer starts (at M2 ~ 21 Mq) before the 
secondary has exhausted hydrogen in the core, and contact is 
avoided (see text). 

4.1.2. Case A 

Fig. ^ demonstrates a completely different mechanism of con- 
tact formation operating in Case A systems, which occurs for 
short periods (cf. also Pols 1994). It shows the radius evolu- 
tion of three secondaries in Case A systems which have identi- 
cal initial parameters except the period (systems No. 53. ..55 in 
Table [3]). All three secondaries evolve through the rapid part of 
the Case A mass transfer (cf. Sect. 3.1) without expanding sig- 
nificantly (which is different for initial mass ratios q < 0.65; cf. 
Sect. 5.1). However, in the two closer systems, the mass trans- 
fer appears so early in the evolution of the primary that the now 
more massive secondary finishes core hydrogen burning first. It 
then starts to expand, but soon fills its Roche lobe and attempts 
a reverse Case B mass transfer (in systems No. 53 and 54). 
Thus, both stars now fill their Roche lobes and come into con- 
tact while the primary is still burning hydrogen in its core; we 
designate this as premature contact. The likely outcome of pre- 
mature contact is a merger because the corresponding periods 
are too small for a successful common envelope ejection. This 
effect marks the lower period limit for contact free binaries. 

In model No. 55, the primary ends core H-burning first. 
Therefore, Case AB mass transfer (cf. Section 3.1) sets in be- 
fore the secondary expands strongly. This has two effects: it 
drives the two stars further apart — the period increases from 
5.24 d at the beginning to 93 d at the end of Case AB mass 
transfer — , and it enhances the mixing of hydrogen into the 
core (cf. Braun & Langer 1995) and thus delays the evolution of 
the secondary. I.e., if the initial period exceeds a critical value. 
Case AB mass transfer starts before the secondary expands and 
premature contact can be avoided. If the secondary star does 



not rejuvenate, reverse Case B mass transfer is avoided as well. 
Curiously, the secondary star in our Case A system No. 1 fin- 
ishes core hydrogen burning just after the end of the Case A 
mass transfer and before the beginning of the Case AB mass 
transfer. However, the accretion due to the Case AB mass over- 
flow leads to a shrinkage of its radius, and contact is avoided 
(see also Fig. 1 1 below). 

The lower period limit for contact free binaries depends 
sensitively on the convection criterion. As mentioned in 
Section 2, we compute stellar models using the Ledoux crite- 
rion for convection and semiconvection with an efficiency pa- 
rameter of Qfsc = 0.01. The semiconvective mixing efficiency 
— which is infinite in the case of the Schwarzschild criterion 
for convection — determines to a large extent whether an ac- 
creting core hydrogen burning star rejuvenates or not (Braun & 
Langer 1995). If it fails, which is the more likely for less ef- 
ficient semiconvection and the more advanced the star is in its 
evolution when it starts to accrete, it retains an unusual struc- 
ture, with a smaller helium core mass and a larger hydrogen 
envelope mass compared to single stars. As such a structure 
keeps the star relatively compact during all the advanced burn- 
ing stages — i.e., it avoids the red supergiant stage (Braun & 
Langer 1995) — it will not fill its Roche lobe after core hy- 
drogen burning. This is the situation in all our Case A systems 
which avoid premature contact and have initial mass ratios of 
g > 0.65. 

For models computed with the Schwarzschild criterion, the 
secondaries in Case A systems rejuvenate during or shortly af- 
ter the rapid Case A mass transfer. This results in longer main 
sequence life times for the secondaries, and thus leads to less 
premature contact situations. I.e., the critical initial period for 
premature contact is shifted to smaller values. In this case, the 
secondaries all expand to red supergiants. Assuming the pri- 
mary star is still present at that time, this is likely to initiate 
reverse Case B mass transfer As the life time of a rejuvenating 
star is extended, the primaries have often finished their evolu- 
tion by that time and the binary may even be disrupted (cf. Pols 
1994). However, as the life time of a low mass helium star (i.e. 
the primary of such a system) can be comparable to the remain- 
ing hydrogen burning life time of the very massive secondary, 
reverse Case B mass transfer onto the low mass helium star 
may also occur. In this case, the helium star builds up a hydro- 
gen envelope and swells to red giant dimensions. We designate 
the ensuing contact situation in Table 3 as reverse contact. 

This scenario is illustrated by five of our Case A systems 
(No. 25, 48, 50, 59, and 72) which are computed using the 
Schwarzschild instead of the Ledoux criterion (marked "c" in 
Table 3). No. 25 and 59, initiaUy having a 12 Mq and a 16 
primary, a mass ratio of 0.67 and 0.69, and a period of 2.5 d 
and 2.0 d, respectively, avoid premature contact. However, they 
encounter reverse contact when the primaries have burnt 9 1 % 
and 80% of their helium in the core. Both systems avoid con- 
tact altogether when they are computed with our standard con- 
vection physics, and evolve a reverse supernova order I.e., an 
increase of the semiconvective mixing efficiency may reduce 
the number of systems with premature contact but increase the 
number with reverse contact. From system No. 25 we conclude 
that, were the Schwarzschild criterion correct, the majority of 
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Case A binaries with 12 Mq primaries would evolve into re- 
verse contact. Comparing our models No. 58 and 59, which 
have identical initial conditions but vary in the treatment of 
convection, shows that also in systems with 16 M0 primaries 
the lower period limit for contact free systems increases. 

In the three other systems computed with the 
Schwarzschild criterion, systems No. 48, 50, and 72, the 
primaries explode as supernovae before the reverse mass 
transfer can occur. Two of these (No. 48 and 72), when 
computed with our standard convection physics, have a reverse 
supernova order (see systems No. 47 and 71 in Table 3), the 
other (No. 50) does not (cf. system No. 49 in Table 3). 

Thus, the lower period limit for contact free evolution with 
the Schwarzschild criterion is expected between the lower pe- 
riod limit and the limit for reverse supernova order for Ledoux 
criterion and Osc — 0.01. The one early case A system con- 
taining a 25 M0 primary (No. 72) also shows no reverse mass 
transfer before the supernova explosion of the primary, despite 
being a system with reverse supernova order when computed 
with our standard convection physics. We conclude that for 
large masses, the amount of contact free case A systems does 
not change very much for different convection physics. On the 
contrary, at 12 primary mass and below no contact free 
case A system are expected at all if the Schwarzschild criterion 
is used. 

4.2. Contact formation due to extreme mass ratios 
4.2.1. Case B 

While the previous two subsections concerned contact forma- 
tion due to large (delayed contact) or small (premature and re- 
verse contact) initial periods, we investigate here contact for- 
mation due to large or small initial mass ratios. 

Fig. ^ shows the mass transfer rate as function of the trans- 
ferred amount of mass for three Case B systems with iden- 
tical initial primary mass (12 M©) and period (6 d), but dif- 
ferent initial secondary masses. For smaller initial mass ratios 
Qi ~ M2;i/Mi^i, the mass transfer rates are somewhat larger 
in the early part of the mass transfer, since the orbit shrinks 
more rapidly for smaller gi. For conservative evolution (also 
ignoring stellar winds), the minimum separation depends on 

initial separation dj and mass ratio as dmin — di ^ (^g^^'iyi ^ ■ 
As, in a given system, period and orbital separation scale as 
(X (P, this means that in systems with the same initial pe- 
riod, the primary star is, at the time when P = Pmin, squeezed 
into a smaller volume for smaller gi. This makes the mass trans- 
fer rate during the first phase of Case B mass transfer larger for 
smaller qi. 

However, the maximum mass transfer rate, which is 
achieved after mass ratio reversal (cf. Sect. 4.1), is very sim- 
ilar in all three systems, even though the post-mass transfer pe- 
riods (115 d, 87 d, and 72 d for secondary initial masses of 
1 1 M0, 9 M0, and 8 Mq, respectively) are not the same. Also 
the amount of mass which is transfeiTed is independent of the 
secondary mass — it consists of the whole hydrogen-rich en- 
velope of the primary. 



. . . , 
: P; = 6 d 


' ' ' \ 


■ ■ , ■ _ 
■ 


; = 


8 Ms (No. 27) 




- M,j = 


9 Ms (No. 19) 




= 


1 1 Ms (No. 3) 





2 4 6 8 10 

"irons [Ms] 



Fig. 8. Mass transfer rates as function of the transfered amount 
of mass for three Case B systems with 12 Mq primaries. The 
initial masses of the secondaries are 1 1 Mq (system No. 3), 
9Mq (No. 19), and 8 Mq (No. 27), as indicated. The initial 
period is 6 d for all three systems. 



Whether or not contact is reached is here mostly determined 
by the reaction of the secondary star to the accretion. Fig. ^ 
shows that the primary with the lowest initial mass swells most 
during the mass transfer. The smallest initial mass ratio for 
which contact is avoided corresponds roughly to the condition 
that the mass accretion time scale of the secondary A/2/M2 re- 
mains smaller than its thermal time scale GM|/(i?2i2)- We 
designate contact due to this reason as q-contact. 

As q-contact is established during the first phase of 
Case B mass transfer, it is quite distinct from delayed con- 
tact (Sect. 4.1). As q-contact evolves early during Case B mass 
transfer, it does not allow a major amount of mass to be ac- 
creted by the secondary. I.e., either the major part of the pri- 
mary's hydrogen-rich envelope can be ejected from the sys- 
tem, or both stars merge. The situation that the secondary star 
can accrete significant amounts after a non-conservative early 
q-contact phase is unlikely since it would bring both stars even 
closer together. 

In principle, there is also a maximum initial mass ratio for 
contact-free evolution of Case B systems: for qi = 1, both stars 
age at the same rate and attempt to expand to red giants simul- 
taneously. However, since the limiting value is very close to 
qi = 1, we ignore this effect here as statistically insignificant. 

4.2.2. Case A 

The maximum mass transfer rate in Case A systems is achieved 
during the rapid mass transfer phase while the orbit is shrink- 
ing. Fig. |l^ shows the mass transfer rates of five Case A sys- 
tems with identical initial conditions except for the initial sec- 
ondary mass and thus the initial mass ratio q-^. As shown in 
Sect. 4.2.1, the smaller qi the larger is the maximum mass trans- 
fer rate since the primary is squeezed into a smaller volume. 
Also, more mass is transfered during the rapid Case A mass 
transfer for smaller q-^. As for Case B, q-contact may occur 



12 



S. Wellstein et al.: Formation of contact in massive close binaries 




P; = 6 d 

— Mji = 8 Mj, (No. 27) 
J = 9 Mg, (No. 19) 
J = 11 Ms, (No. 3) 



16 
14 

12 



^ 10 



2.5 d 

1 1.5 Mg (No. 1) 
10.5 Us (No. 10) 
■ M,; = 9.5 Mg (No. 15) 
M,; = 8.5 Uj, (No. 23) 
= 7.5 Me (No. 31) 



20 



22 
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Fig. 10. Mass transfer rate as a function of transfered amount 
of mass, for five Case A systems with the same initial primary 
mass (12 Mq) and initial period (2.5 d). The initial secondary 
masses are 11.5 Mq (System No. 1), 10.5 M© (No. 10), 9.5 Mq 
(No. 15), 8.5 Mq (No. 23), and 7.5 Mq (No. 31). 



for small qi due to higher maximum mass transfer rates and 
due to larger thermal time scales of the secondary, the latter 
effect being the more important one. In Fig. |ll], all systems 
shown avoid q-contact, but the relative radius increase during 
the rapid Case A mass transfer is clearly larger for the systems 
with smaller initial secondary masses. 

The second expansion of the secondaries during the Case A 
mass transfer (Fig. [ll]) is due to their nuclear evolution dur- 
ing the slow mass transfer phase. It is ended by the onset of 
Case AB mass transfer. The secondary star in system No. 1 
has already ended core hydrogen burning and expanded appre- 
ciably before the onset of Case AB mass transfer The onset of 
Case AB mass transfer stops its expansion and prevents the sys- 
tem from evolving into premature contact. However, Figs. |lo| 
and [n] demonstrate, that for too large mass ratios, massive 
Case A binaries evolve into premature contact as in the case of 
too small initial periods (cf. Section 4.1.2). I.e., while Case B 



Fig. 11. Radii of the secondary stars as function of their mass 
for the same systems shown in Fig. |l0[ The first maximum is 
achieved during rapid Case A mass transfer. The second max- 
imum is achieved during slow Case A mass transfer due to the 
nuclear evolution of the secondary. The radius decrease after 
the second maximum is due to the Case AB mass transfer. 



systems evolve into contact only when their initial mass ratio is 
too large. Case A binaries can do so for too large or too small 
initial mass ratio. 
Fig. 
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shows also, that the total amount of mass which is 
transfered during case A and AB is nearly independent of q^. 
This can not be expected in general, as we shall see below that 
the primary mass after Case A and AB mass transfer depends 
on the core hydrogen abundance at the onset of mass transfer 
and on the initial mass ratio (cf. Fig. below). While smaller 
initial mass ratios lead to smaller primary masses, the larger 
Roche lobe due to the lower secondary mass (for the same pe- 
riod) leads to mass transfer only later in the evolution of the 
primary and thus to larger post-mass transfer primary masses. 
Both effects cancel each other in the models shown in Fig. |l^. 

5. Discussion 

5. 1. The contact-free regime 

Fig. |l2| and |r3| summarise the range of contact-free evolution 
for systems with 12 Mq and 16 Mq primaries, respectively. For 
all mass transfer cases and primary masses, there is a rather 
well defined critical initial mass ratio qc such that only sys- 
tems with larger q^ avoid contact. For q^ < q^, they evolve into 
q-contact. Likely, most of the systems undergoing g-contact 
merge, since contact occurs early during the mass transfer, at 
a small orbital separation. We find q^ ~ 0.65 for all Case A 
systems with primaries in the initial mass range 12...25Mq. 
For Case B systems, qc decreases from qc ~ 0.7 at 12 Mq to 
qc ~ 0.8 at 16 Mq and qc ~ 1 al 25 Mq initial primary mass. 

The range of contact-free Case A systems is further limited 
by the possibility of premature contact, which occurs for ini- 
tial periods below a certain threshold, but also for initial mass 
ratios above a critical value, i.e. too close to one. We should 
note that for even smaller initial periods than investigated here. 
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premature contact may occur already before the secondary has 
finished core hydrogen burning (cf. Pols 1994). In any case, 
the outcome of premature contact is most likely the merging of 
both stars. The possibility of reverse contact and of a reverse 
supernova order for the Case A systems, and their dependence 
on the semiconvective mixing time scale are discussed in the 
next section. 

As described in Section 4. 1 .2, the uncertain efficiency of 
semiconvective mixing introduces a considerable uncertainty 
to the evolution of Case A systems. Semiconvective mix- 
ing affects the rejuvenation process in an essential way, as 
it is this process which transports hydrogen into the convec- 
tive core of the accretion stars. While for infinitely fast semi- 
convective mixing — which corresponds to the use of the 
Schwarzschild criterion for convection — all accretion stars 
rejuvenate (Hellings 1983), this is not the case for finite semi- 
convective mixing time scales (Braun & Langer 1995). Longer 
mixing time scales lead to shorter core hydrogen burning life 
times of the secondaries. Therefore, slow semiconvection — 
as assumed here — leads to earlier contact. Nevertheless, slow 
semiconvection leads perhaps not to a smaller total number of 
contact systems, since for those systems which avoid prema- 
ture contact, the chance of reverse contact during core helium 
burning of the secondary is much smaller than in the case of 
fast semiconvective mixing (cf. Sect. 4.1.2). 

We note that although semiconvection can be important 
for the upper end of the mass spectrum of our study (cf. also 
Wellstein & Langer 1999), another effect dominates for the less 
massive Case A systems. In our systems with an initial primary 
mass of 12 Mq, the luminosity ratio after the mass transfer is 
very large, of the order of 100, while it is only of the order 
of 10 for our most massive systems (see Table 3). The reason 
for the drop of the luminosity ratio for higher primary masses is 
simply the flattening of the stellar mass-luminosity relation for 
increasing mass. As the luminosity ratio reflects the ratio of the 
evolutionary speeds of both components, it is clear that with a 
post-mass transfer value of 100, only very late Case A systems 
can avoid contact — independent of whether the secondary re- 
juvenates or not. The large luminosity ratio for lower primary 
initial masses is also the main reason why the initial parame- 
ter space for contact-free evolution becomes smaller for lower 
mass systems. 

Finally, there is an upper period which limits the contact- 
free regime of the Case B systems. Our calculations have 
shown that systems close to this limit may experience delayed 
contact, with the consequence of considerable accretion on the 
secondary, and a high chance to avoid the merging of both stars. 
Initially even wider systems may correspond more to the stan- 
dard common envelope scenario, again with the possibility to 
avoid a merger, but likely without significant accretion of the 
secondary. 

By comparing Figs. |l2| and [l^, and including our results 
for even more massive systems (see Table 3, and Wellstein 
& Langer 1999), it becomes clear that the range of contact- 
free Case A systems becomes smaller for smaller initial pri- 
mary masses, while that for contact-free Case B systems be- 
comes smaller for larger ones. E.g., there are no contact-free 
Case B systems for initial primary masses of 25 M©. We con- 
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Fig. 12. Distribution of all computed binaries with 12 M© pri- 
maries in the log P versus q-diagram. Asterisks mark contact- 
free systems, while squares mark systems which evolve into 
contact. Systems marked with triangles are borderline cases, 
i.e., they evolve into a short contact phase but the secondary ra- 
dius never exceeds its Roche radius by more than a factor 1.5. 
The soUd line separates Case A (below) and Case B systems. 
All case A systems for this primary mass have a reverse su- 
pernova order The dashed lines indicate the boundary between 
contact-free and contact evolution. The dashed-dotted line is 
defined by the condition that the primary fills its Roche lobe 
already on the zero age main sequence. 



1.2 






1 1 1 1 




1 ' ' 'd 


1.0 








□ 


A _ _ - - 


0.8 




□ 


□ 


I 

1 ^ 




0.6 




□ 


A 








A 1 %. 














* ... 






0.4 
0.2 


1 ^ 

u - ' 


□ 
□ 




□ 


^ - -nr " " - 


0.0 


1 1 1 1 



0.50 



0.60 



0.70 0.80 

q 



0.90 



1.00 



Fig. 13. The same as in Fig. [ij but for systems with initial pri- 
mary masses of 16 Mq. The dotted line divides the contact-free 
Case A systems into such with normal and reverse supernova 
order, the latter occurring in the systems below the line. 



elude that at lower masses, conservative evolution is to be 
expected mainly from close Case B systems, while at larger 
masses Case A systems dominate the conservative evolution. 
Considering that massive binaries with initial periods of up to 
some thousand days may interact, and adopting an equal num- 
ber of binaries per log P bin, up to one third of all interacting 
systems in the considered mass range and with qj > 0.7 may 
evolve conservatively. 
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Fig. 14. JVIasses after Case AB mass transfer for primaries with 
an initial mass of 16 M©, as function of their central helium 
abundance at the onset of Case A mass transfer For compari- 
son, the post Case B mass transfer masses of 16 Mq primaries 
are all about 3.6 M©. 

5.2. Remnant types 

Those of our systems which avoid merging evolve into 
OB+He star systems. We discuss in the following the further 
evolution of the primaries, the secondaries, and of merger stars. 

5.2.1 . Primaries 

Table 3 shows that the helium star masses after the Case AB/B 
mass transfer in our systems fall into the range 1...7 M©. The 
helium star masses in the Case B systems correspond simply to 
the helium core mass of the primary at core hydrogen exhaus- 
tion, with a very small scatter (+0.1 M©). Fig. |l^ shows that 
the situation is more complex for the Case A systems. We find a 
strong dependence of the helium star mass on the initial period, 
i.e. on the central hydrogen abundance of the primary when the 
mass transfer begins. We also find a weaker dependence on the 
initial mass ratio (cf., systems No. 41, 55, 58). Fig. ^demon- 
strates that our 16 M© Case A primaries become 2. ..3 M© he- 
lium stars — compared to a Case B post-mass transfer mass of 
~ 3.6 M©, which is independent of other parameters. We point 
out that these strongly reduced Case A primary masses are of- 
ten ignored in simplified binary studies which rely on single 
star evolutionary models. 

The further evolution of the helium stars depends strongly 
on their mass, and partly on the evolution of their companion 
stars. The more massive helium stars in our mass range corre- 
spond to Wolf-Rayet stars, and will undergo wind mass loss be- 
fore they explode as Type Ib/Ic supernova. Those helium stars 
with masses below ^ 4 M© expand to giant dimensions after 
core helium exhaustion. Depending on whether the companion 
star is still present at that time, this may lead to Case ABB/BB 
mass transfer, and a further reduction of the mass of the he- 
lium star by ~ 0.5 M© or so (cf. system No. 37 in Sect. 3.2). 
Independent of this, the helium stars with final masses in the 
range 2...4M© are likely to explode as helium giants. The ef- 



fect of the large radii of these stars on the supernova light curve 
is still unexplored. Helium stars in the mass range 1.4...2M© 
could, in principle, also make it to core collapse, but as their 
post-core helium burning life time is rather large it may be 
more likely that stellar winds will reduce their mass below the 
Chandrasekhar mass and they end up as white dwarfs, like their 
counterparts with masses between 1...1.4M©. 

We point out that all 12 M© primaries in our contact-free 
Case A systems end up as white dwarfs, and even some 16 M© 
primaries may do so (e.g., system No. 55). I.e., the limiting 
initial mass separating white dwarf formation from core col- 
lapse is largely increased in Case A binaries, compared to sin- 
gle stars. We note that this effect exists also for Case B binaries, 
strongly so for primaries which undergo a Case BB mass trans- 
fer, but at a moderate level even otherwise. The same effect for 
the limiting initial mass between neutron star and black hole 
formation has been discussed by Wellstein & Langer (1999). 

Finally, we can not exclude that some of the helium stars 
with masses close to but above the Chandrasekhar mass de- 
velop a degenerate C/O-core and ignite carbon explosively. In 
that case, they would produce a Type la supernova. However, 
in the 1.45 M© helium star of system No. 31 (cf. Sect. 3.1), a 
luminosity of ^ 10"* L© during shell helium burning implies a 
growth rate of the C/O-core mass of - 10"^ M© yr~\ which 
is fast enough to allow carbon to ignite non-explosively. 

5.2.2. Secondaries 

All our secondary stars evolve to core collapse. In Case A sys- 
tems, about half of the secondaries does so before the primary 
has ended its evolution (cf. Table 3). As most of our primaries 
did not rejuvenate, they remain blue supergiants throughout 
their post-main sequence evolution (cf. Figs. 1 and 3). When 
they explode as Type II supernova, they will thus resemble 
SN 1987A. We noted in Sect. 3.1 that the secondary of sys- 
tem No. 31 provides a viable progenitor model for SN 1987A. 
In this case, there should be a 1.5 M© helium star left in the 
supernova remnant. In similar systems, however, the supernova 
order is not reversed and no such helium star would be expected 
(cf. Table 3). 

5.2.3. IVIerger 

Several of our systems evolve into contact and are likely to 
lead to a merger While it is beyond our scope to predict the 
properties of the merged stars, we can distinguish two types of 
events. If the merging occurs during Case A mass transfer and 
before the secondary has terminated core hydrogen burning, 
the merged object will still be a core hydrogen burning star. 
Its surface may be enriched in hydrogen burning products, but 
otherwise it might resemble a normal main sequence star in 
many respects. 

However, a merger during a Case B or reverse Case B mass 
transfer (cf. Sect. 4.1.1; note that the latter may even occur dur- 
ing a Case A mass transfer; cf. Fig. 7) has a different result. 
Since the original primary has already developed a compact 
hydrogen-free core, the merger star will have a hydrogen-free 
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core with a mass close to that of the primary. It is uncertain how 
much of the hydrogen-rich material is lost during the merging 
process, and whether some hydrogen-rich matter can penetrate 
into the helium core. However, if the mass of the merger star 
exceeds the mass of the primary star, which appears likely, then 
the helium core mass of the merger star is much smaller than 
the helium core mass of a single star of comparable mass. I.e., 
its structure is similar to that of our non-rejuvenating accretion 
stars during core helium burning. As a consequence of their 
small helium core mass, these stars may avoid to evolve into 
red supergiants (Podsiadlowski et al. 1990, 1992). 

5.2.4. Reverse supernova order 

In our Case A binaries, a reversal of the supernova order ap- 
pears in all systems containing a 12 Mq primary (c.f. Fig.Jo) 
and in some systems containing a 16 Mq primary (cf. Fig. pj^ . 
At 12 Mq, supernova order reversal depends not so much on 
whether or not the secondary star rejuvenates. Instead, the mass 
ratio after the rapid phase of the Case A transfer, and thereby 
the luminosity ratio, is essential. Due to the strong decrease of 
the primaries mass, it becomes a factor of about 100 less lumi- 
nous than the secondary. Therefore, independent of how much 
fresh hydrogen is mixed into the secondaries core, it finishes its 
evolution first as long as no reverse mass transfer occurs. 

However, in case of rapid semiconvective mixing, the sec- 
ondary may rejuvenate, with the consequence of reverse mass 
transfer before it can evolve into a supernova. Due to the ex- 
treme mass ratio, the reverse mass transfer is expected to be 
non-conservative. Its outcome is either a merger or a short pe- 
riod double helium star binary. Assuming the latter, we can es- 
timate whether or not the supernova order would be reversed. In 
system No. 17 — computed with age — 0.02 — the secondary 
is more evolved at time of reverse mass transfer, i.e. both stars 
burn helium in their cores, with central helium mass fractions 
of 0.5/0.3 for the primary/secondary. Thus the supernova order 
is expected to be reversed. For system No. 25 computed with 
Schwarzschild criterion, it is less clear what happens. The pri- 
mary has already spent about 6 10^ yr at central helium burning 
and has reached a core helium mass fraction of 0.09. The sec- 
ondary has reached a core helium mass fraction of 0.63 after 
only 3 10^ years of central helium burning. Because the pri- 
mary's helium shell burning phase is expected to last more than 
1 10^ years (cf. system No. 31 in Table [l|) it appears possible 
that also in this system the supernova order is reversed, if the 
system does not merge. 

For systems with an initial primary mass of 25 Mq the 
convection physic dominates the question whether the super- 
nova order is reversed. This can be seen in systems No. 65, 66 
and 67, where the value for the efficiency parameter for semi- 
convective mixing age is the only difference. Only in system 
No. 65 with asc — 0.01 the supernova oder is reversed. With 
ttsc > 0.02 no reversal appears. Also the case B system 68 re- 
verses its supernova order with a^c = 0.01. This shows that 
the mass of the primary after mass transfer in this systems is 
not very important because the nuclear time scales of the stars 
in this mass range vary less with mass than for lower mass 



stars. Also the early case A binary No. 71 has a reverse su- 
pernova order in contrast to system No. 72, which has the same 
parameters, but is calculated with the Schwarzschild criterion. 
The systems with 16 Mq primary are an intermediate case, but 
seems to behave more like the systems containing a 25 Mq 
primary (c.f. systems No. 47 and 48). 

Whether or not a reverse mass transfer occurs before the 
primary becomes a compact object depends strongly on the 
semiconvective mixing speed. For fast semiconvective mixing, 
the secondaries become red supergiants. Then, the crucial ques- 
tion is not if reverse mass transfer occurs but when it occurs. A 
shorter semiconvective mixing time scale has two effects. First, 
the maximum radius of the secondary after core hydrogen ex- 
haustion increases (cf. Wellstein & Langer 1999). This makes 
it more likely for the secondary to fill its Roche lobe and trans- 
fer matter back to the primary. Second, the time the secondary 
spends on the main sequence increases, as more hydrogen is 
mixed into the convective core during central hydrogen burn- 
ing. Thus it is less likely that the reverse mass transfer occurs 
before the primary has finished its evolution and has become a 
compact object. 

Although in general, systems which have a normal super- 
nova order would be able to avoid reverse mass transfer and 
vice versa, some systems which show a reverse supernova or- 
der for slow semiconvective mixing (ctsc — 0.01) avoid re- 
verse mass transfer during the lifetime of the primary when the 
Schwarzschild criterion is used, as shown by systems No. 48 
and 72. 

5.3. Comparison to previous work 

Pols (1994) presented 24 binary models for initial primary 
masses in the range 8...16Mq, mostly for Case A systems. 
Since he used the Schwarzschild criterion for convection, the 
secondary components rejuvenated in all computed cases. As a 
consequence, contact phases during slow case A mass transfer 
occurred for shorter periods than in our models. Furthermore, 
due to the rejuvenation, a supernova explosion of the secondary 
as the first star could only appear after a reverse mass transfer 
phase. This mass transfer appears at very large mass ratios and 
must therefore be non-conservative. As discussed by Pols, it 
is an open question whether or not such a binary merges or is 
able to expel the secondary's hydrogen-rich envelope to form a 
close binary consisting of two helium stars. As the helium star 
resulting from the secondary is then the more massive compo- 
nent, a supernova order reversal is possible if the merging of 
both stars can be avoided. 

Sybesma (1985) computed 5 binary sequences starting with 
a primary mass of 20 Mq, a secondary mass of 10 Mq, and 
periods ranging from 1.5 up to 10 days. His results differ 
from ours mainly due his assumption of huge convective core 
overshooting. This leads to much more massive helium-rich 
cores in the main sequence models. Consequently, he obtains 
Case A mass transfer for initial periods which fall in our Case B 
regime. As Case A mass transfer is more likely conservative in 
this mass range, Sybesma obtains conservative evolution for 
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larger periods than we (e.g. 10 days), and for smaller initial 
mass ratio (q=0.5). 

Several authors specify the fraction (3 of the primary's enve- 
lope which is accreted by the secondary in advance. E.g., based 
on observations, De Greve & de Loore (1992) and de Loore 
& De Greve (1992) use /3 = 0.5, which means half of the 
primary's envelope mass leaves the system, together with an 
amount of angular momentum which needs to be specified ad- 
ditionally. While an average value of /? = 0.5 appears in fact 
possible, our models favour the idea that (3 decreases continu- 
ously for larger initial periods, with (3=1 for systems which 
avoid contact and (3 = for large initial periods. E.g., we ob- 
tained (3 = 0.68 and (3 = 0.60 for our systems No. 7 and 8 
(cf. Table 3). By comparison with De Greve & de Loore (1992) 
and de Loore & De Greve (1992), we note that the contact-free 
regime is expanded somewhat if (3 = 0.5 is used ab initio. 

6. Comparison with observations 

While it is unlikely that any of the binary systems considered 
here has an observable counterpart during the mass transfer 
— with the exception of the slow phase of the Case A mass 
transfer — they have a relatively long-lived post-mass transfer 
stage. We emphasise that this is true for those binary which sur- 
vive the mass transfer, but also for those which lead to mergers. 
Therefore, we characterise the typical properties of systems af- 
ter the first mass transfer phase in Table 4. At this stage, the sys- 
tems consist of a main sequence or supergiant O or early B star 
of 18. ..40 M0, and a helium star in the mass range 1...7 Mq. 
The orbital periods range from 24 d to 255 d. 

An inspection of Table 4 indicates that it may be a chal- 
lenging task to try to detect systems like these observationally. 
A photometric investigation will have limited success due to 
three factors. First, the luminosity ratio of both stars is large, 
mostly of the order of 10, with the worst case of 160 (system 
No. 24) and the best of 4 (system No. 68). Second, the less 
luminous star has a much smaller radius than the brighter star, 
which minimises the brightness amplitude in case of an eclipse. 
Third, both stars underfill their Roche volumes by a large fac- 
tor, which means that they are Ukely to radiate isotropically, 
and that eclipses are not likely. 

We note that after core helium burning, all the helium 
stars below ~ 4 Mq expand to giant dimensions and — if the 
OB star has not yet turned into a supernova — may even fill 
their Roche volumes (cf. Sect. 3.2). As also their luminosities 
increase by up to a factor 10, a photometric detection may in- 
deed be more likely in this stage. However, for the more mas- 
sive helium stars (M[ J> 3 M0), this stage is very short lived 
(^ 0.1 Myr) compared to their core heUum burning Ufe times 
(~ 1 Myr; cf. Table 2). In systems with lower mass helium 
stars, the time scales are more favourable: the 1.45 Mq helium 
star of system No. 3 1 burns helium in its core for about 5 Myr, 
and its ensuing shell burning stage lasts at least 1 Myr, possi- 
bly more depending on the mass loss rate during this stage (cf. 
Sect. 3.1 and Table 1). 

Table 4 also shows the orbital velocities of both com- 
ponents after the mass transfer. While the helium stars have 
favourably large velocities of the order of 100. ..200 km s~^. 



they might be unobservable against the glare of their OB com- 
panion. The latter moves with only 5. ..35 km s~^. While these 
velocities are small, it may not be hopeless to try to detect these 
motions in OB stars. 

It is interesting to note that the assumption of conservative 
evolution leads, in all respects, to system properties which are 
less favourable for an observational detection of such systems, 
compared to the assumption of mass- and angular momentum 
loss from the system. This is demonstrated by the estimated 
outcome of our delayed contact system No. 7 in Sect. 4.1.1, 
which results in a shorter orbital period (7 d), correspondingly 
larger orbital velocities, and a mass- and luminosity ratio closer 
to 1, compared to similar conservative systems. This tendency 
can also be verified in the results obtained by de Loore & 
De Greve (1992) obtained with (3 = 0.5, i.e. assuming that 
half of the transferred mass leaves the system. 

The best chances to find post-mass transfer systems like 
those of Table 4 seems to be at the upper end of the consid- 
ered mass range. Our systems with initial primary masses of 
25 M0 obtain the smallest luminosity ratio L'2fL[, the short- 
est periods, and the largest orbital velocities. The helium star 
masses in these systems are in the range 5. ..7 Mq, which makes 
the stars comparable to Wolf-Rayet stars, or, more specifically, 
WN stars. However, although we know 12 OBh-WN systems 
with established mass ratios and orbital periods (van der Hucht 
2000), none of them fits to our models. While 8 out of these 
12 binaries have clearly much higher initial system masses 
than those systems considered here, the remaining four (WR3, 
WR31, WR97, and WR139; cf. van der Hucht 2000) all have 
mass ratios of Mqb/Mwr ^ 3, and mostly much shorter peri- 
ods than our models. 

There are two possible ways to interpret the disagreement 
between our models and the observed OBh-WN binaries. One 
possibiUty is that our assumption of conservative evolution 
does not hold for systems with primaries of initially ~ 25 Mq. 
This can be motivated by the increasing role of radiation pres- 
sure in higher mass systems (cf., Schuerman 1972; Vanbeveren 
1978). On the other hand, strong non-conservative effects in 
massive systems may result in an uncomfortably large upper 
critical ZAMS mass for neutron star formation of ^ 40 Mq, as 
imposed by the massive X-ray binary Wray 977 (Kaper et al. 
1995, Ergma & van den Heuvel 1998, Vanbeveren 1998b). 

The other possibility involves to estimate the probabil- 
ity of finding one system of the kind computed here among 
four OBh-WR binaries. As outlined in Sect. 5.1, the param- 
eter range for conservative evolution becomes narrower for 
higher masses. For our systems with 25 Mq primaries, contact 
is avoided for an initial period of 4 d, while 5 d resulted al- 
ready in contact (cf. systems No. 65. ..74 in Table 3).Assuming 
initially an equal number of binaries per log P-interval, and an 
upper period limit for interaction of 2000 d, results in a ratio of 
contact- free versus contact surviving binary systems of 1/3.7 if 
all Systems with initial mass ratios q < 0.6 merge. If only the 
systems with initial mass ratios q < 0.6 and initial periods P 
< 200 d merge this ratio becomes 1:6.6. Because it is not well 
known how much systems really merge — also an unknown 
fraction of the large period contact systems might actually lead 
to a merger rather than to a OBh-WR binary — we assume that 
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Table 4. Observable quantities for our contact-free systems. After system number, initial values of primary mass Mi, secondary 
mass M2, and period Pj, we give quantities obtained after the Case AB or Case B mass transfer, at a time when the core helium 
mass fraction of the primary has decreased to 0.8 during core helium burning: primary and secondary mass and mass ratio, 
bolometric luminosities, effective temperatures (where no wind effects are included for the helium stars), stellar radii, the radius 
of the secondary in units of its Roche radius, orbital separation, system period, and orbital velocities. As in Tab. 3, bold values 
for the initial masses identify the component which ends its evolution first, typically by a supernova explosion. Both components 
of system No. 43 explode simultaneously, within the accuracy of the computation. 
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the ratio is between 1:4 and 1:7. Thus the probability of finding 
one suitable binary among four may be of the order of 1/2 or 
lower. However, the survivors of the non-conservative contact 
evolution end up as shorter period systems with a mass and lu- 
minosity ratio closer to one[]than the contact-free systems, and 
are thus detected much easier. 

In summary, the failure to detect counterparts of our more 
massive post-mass transfer systems is inconclusive. For our 
lower mass systems, i.e. such with initial primary masses of 

^ This fits well to the properties of the WR binaiies WR31, WR97, 
and WR139, while WR3 has most likely not experienced any mass 
transfer 



12 M0 and 16 Mq, the chance to detect counterparts may even 
be lower than for our higher mass systems (cf. Table 4). In fact, 
we do not know of any suitable OB star with a 1...3.6 M0 he- 
lium star companion. The object which comes closest to our 
models may be $ Persei, a ^ 9 B star with a 53 000 K hot 
1.1 M0 companion in a 127 d orbit (Gies et al. 1998). Its prop- 
erties, and the fact that the B star is spinning rapidly, support 
actually the possibility of a close-to-conservative evolution at 
least in this system. A more puzzling system which may be 
relevant in the context is A Sco, a ^ 10.5 Mq B dwarf with a 
massive white dwarf companion in a ~ 6 d orbit (Berghofer et 
al. 2000). Waters et al. (1989) speculate that unseen compan- 
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ions in a number of Be stars may actually be helium stars, in 
accordance with Pols et al. (1991) and Portegies Zwart (1995). 

The lack of OB+He star systems which are more massive 
than $Per but less massive than OB+WR systems is, any- 
how, very remarkable. I.e., if the evolution of binaries with 
12. ..16 M0 primaries would, for a big fraction of the possible 
period range, proceed non-conservatively and without leading 
to mergers, we would expect short period OB+He star systems 
which are easier detectable, with a smaller mass- and luminos- 
ity ratio compared to our models. Therefore, a non-detection 
of such systems may be easiest understood if the conservative 
evolution makes up for a large fraction of the non-merging in- 
teracting binaries in the considered mass range. 

Whichever way it is, there are two kinds of observations 
hinting for a large number of OB+He star binaries waiting to 
be discovered. One is the existence of massive X-ray bina- 
ries, which are likely evolving from an OB+He star stage (e.g., 
Moers & van den Heuvel 1989). The second is the relatively 
large number of Type Ib/Ic supemovae, which are mostly at- 
tributed to exploding helium stars with masses above 2M0 
(Podsiadlowski et al. 1992). I.e., it might be a rewarding project 
to attempt an observational search for helium star components 
to OB stars. In that respect. Table 4 can give indications of the 
expected parameters, with the remark that if mass transfer is 
partly non-conservative, then mass ratio, luminosity ratio, pe- 
riod and orbital velocities become more favourable for a detec- 
tion. 

Assuming a ratio of interacting binaries to single stars of 
0.3, and ignoring the unknown fraction of mergers among the 
binaries, implies that 1 out of 3 OB stars in the considered mass 
range would, during its hfe time, obtain a He star companion. 
Depending of the rejuvenation of the secondary, this phase may 
last a considerable fraction, say one third, of the secondary's 
life time. This optimistic estimate indicates that only one out of 
ten OB stars in the considered mass range may have a helium 
star companion. However, these OB stars have a surface com- 
position which is typically enriched in hehum by 10. ..20%, and 
which has a N/C-ratio of about one order of magnitude above 
solar (cf. Sect. 3, and Wellstein & Langer 2001). I.e., even 
though also single stars may become hehum- and nitrogen-rich 
during their main sequence evolution due to rotational mixing 
(Heger & Langer 2000), a pre-selection of chemically enriched 
OB stars may lead to a large detection probability of helium 
star companions. 

Another way to pre-select OB stars with potential helium 
star companions is to identify OB type blue stragglers in young 
clusters or associations. The accretion of substantial amounts 
of matter by the primary may lead it to achieve a mass which is 
already above the turn-off mass of the cluster or associations. 
E.g., the 7.5 M0 Case A secondary shown in Fig. 1 grows to 
about 18 Mq at an age of 13 Myr, which is the turn-off age for 
13...14M0. The 15 Mq Case B primary, upon turn-off of its 
16 M0 companion, even grows to 27 M0. The remaining main 
sequence life time of the secondaries after accretion is of the 
order of 1 Myr, which may be long enough to catch them during 
this stage. Furthermore, the accretion of angular momentum 
may spin up the star, and lead to internal mixing processes and 
thereby to an overluminosity of the star with respect to its mass. 



In summary, enriched surface abundances, blue straggler 
nature, overluminosity for the actual mass, and rapid rotation 
in an O or early B star may be hints to a hehum star companion. 

7. Summary and conclusions 

In this paper, we have, on the basis of evolutionary models with 
up-to-date input physics, investigated the regime of contact- 
free evolution in close binaries with initial primary masses in 
the range 12. ..25 M0, under the assumption that contact-free 
systems evolve conservatively, i.e., that no matter except stel- 
lar winds leaves these systems. Our results confirm earlier ob- 
tained ones quahtatively (e.g., Pols 1994), but differ quantita- 
tively and extend over a wider parameter range. 

We demonstrate that the number ratio of contact free 
Case A-to-Case B systems is a strong function of the system 
mass. At the high mass end (25 Mq initial primary mass), 
the contact-free Case B channel is closed completely, while 
contact-free Case A systems exist for a wide parameter range 
(cf. Table 3), although we find this range to depend strongly on 
the semiconvective mixing efficiency (cf. Sect. 5, and Wellstein 
& Langer 1999). In contrast, the contact-free Case A channel 
becomes narrow at 12 Mq initial primary mass, whatever the 
assumptions on semiconvective mixing are, while the contact- 
free Case B range becomes large (Fig. 12). At an initial primary 
mass of 16 M0, the situation is intermediate (Fig. 13). 

We investigate the vastly different types of contact situa- 
tions of systems in the considered mass range, and determine 
the critical periods and mass ratios to avoid contact (Sect. 5.1). 
We also estimate the consequences of the various kinds of con- 
tact. As a new phenomenon, we find that in what we call de- 
layed contact, Case B systems can transfer the majority of the 
primary's envelope conservatively before — towards the end 
of the Case B mass transfer — contact and possibly a common 
envelope evolution occur. We find the delayed contact regime 
may lead, with increasing initial period, to a smooth transition 
from conservative contact-free evolution to the classical com- 
mon envelope evolution, with the parameter (3, which describes 
which fraction of the primary's envelope is accreted by the sec- 
ondary, varying smoothly from /3 = to /3 = 1. We are able 
to compute values of (3 for our delayed contact systems, which 
we find in the range /3 = 0.4. ..0.7. 

How far or close our model sequences actually are to the 
evolution of binary systems in nature depends to a large ex- 
tent on the question whether contact-free systems do indeed 
evolve conservatively. I.e., the accreting star is likely to be 
rapidly spun-up to close to critical rotation (e.g.. Packet 1981), 
and it depends on unknown viscous time scales whether or 
not it will be able to continue to grow in mass (Paczynski 
1991). This question can in principle be answered by com- 
paring our post-mass transfer models with corresponding ob- 
served binaries. However, as outlined in Sect. 6, even though 
virtually no observed counterparts to our post-mass transfer 
models — i.e., no OB+He star systems with large mass ra- 
tios and orbital periods — are known, we could not conclude 
that conservative evolution in massive binaries is not reaUsed 
in nature, for two reasons. First, identifying rather far and faint 
hehum star companions to bright OB stars may be too diffi- 
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cult. Second, non-conservative evolution predicts even easier 
detectable OB+He star systems but, again, no such systems are 
known, at least not corresponding to primary masses between 
12 and 20 Mq — while three known OB+WN binaries might 
indeed be the result of non-conservative binary evolution with 
primaries above 20 Mq. 

We want to emphasise that an answer to the question 
whether massive binaries can evolve conservatively is impor- 
tant to understand the late evolutionary stages: e.g., Type Ib/Ic 
supemovae, massive X-ray and black hole binaries, and 7- 
ray bursts due to massive binaries. In principle one may con- 
strain the conservativity of massive binaries by observations 
of the advanced stages. However, assumptions on even more 
uncertain processes have to be made in such attempts — 
e.g., on common envelope evolution, neutron star birth kicks, 
etc. — which strongly limits the success of such approaches. 
Furthermore, often, e.g. in population synthesis studies, only 
statistical information is obtained. For example, Dalton & 
Sarazin (1995) find the period distribution of massive X-ray 
binaries consistent with an average loss of ~ 70% of the trans- 
ferred envelope mass. However, this can neither support nor 
defeat the existence of conservative evolution for contact-free 
massive binaries. 

The very best way to constrain the first major mass transfer 
phase in massive close binaries would indeed be to identify a 
sample of helium star companions to OB stars. We repeat that 
those are predicted by conservative and by non-conservative 
binary evolution models, but with different periods and mass 
ratios. They are also indirectly proven to exist abundantly by 
the relatively large number of Type Ib/Ic supemovae and mas- 
sive X-ray binaries. In Sect. 6, we identify several ways to pre- 
select OB star candidates with an improved probability for a 
heUum star companion. 
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